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Abstract 

Cloaking devices are prescriptions of electrostatic, optical or elec- 
tromagnetic parameter fields (conductivity a(x), index of refraction 
n(x), or electric permittivity e(x) and magnetic permeability n(x)) 
which are piecewise smooth on M 3 and singular on a hypersurface S, 
and such that objects in the region enclosed by £ are not detectable to 
external observation by waves. Here, we give related constructions of 
invisible tunnels, which allow electromagnetic waves to pass between 
possibly distant points, but with only the ends of the tunnels visible 
to electromagnetic imaging. Effectively, these change the topology of 
space with respect to solutions of Maxwell's equations, corresponding 
to attaching a handlebody to M 3 . The resulting devices thus function 
as electromagnetic wormholes. 
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1 Introduction 



There has recently been considerable interest, both theoretical [HI [HI HH El 
[13] and experimental [20], in invisibility (or "cloaking") from observation by 
electromagnetic (EM) waves. (See also jTTJ for a treatment of cloaking in the 
context of elasticity.) Theoretically, cloaking devices are given by specifying 
the conductivity o~{x) (in the case of electrostatics), the index of refraction 
n(x) (for optics in the absence of polarization, where one uses the Helmholtz 
equation), or the electric permittivity e(x) and magnetic permeability fj,(x) 
(for the full system of Maxwell's equations.) In the constructions to date, 
the EM parameter fields ( a;n;e and fi ) have been piecewise smooth and 
anisotropic. (See, however, [SI Sec. 4] for an example that can be interpreted 
as cloaking with respect to Helmholtz by an isotropic negative index of re- 
fraction material.) Furthermore, the EM parameters have singularities, with 
one or more eigenvalues of the tensors going to zero or infinity as one ap- 
proaches from on or both sides the cloaking surface S, which encloses the 
region within which objects may be hidden from external observation. Such 
constructions might have remained theoretical curiosities, but the advent of 
metamaterials\\\ allows one, within the constraints of current technology, to 
construct media with fairly arbitrary e(x) and n(x). 

It thus becomes an interesting mathematical problem with practical signif- 
icance to understand what other new phenomena of wave propagation can 
be produced by prescribing other arrangements of e and /i. Geometrically, 
cloaking; can be viewed as arising; from a singular transformation ofM 3 . In- 
tuitively, for a spherical cloak [61 [71 [18], it is as if an infinitesimally small 
hole in space has been stretched to a ball D; an object can be inserted in- 
side the hole so created and is then invisible to external observations. On 
the level of the EM parameters, homogeneous, isotropic parameters e, fi are 
pushed forward to become inhomogeneous, anisotropic and singular as one 
approaches S = dD from the exterior. There are then two ways, referred to 
as the single and double coating in [3J, of continuing e, fi to within D so as to 
rigorously obtain invisibility with respect to locally finite energy waves. We 
refer to either process as blowing up a point. As observed in [3J, one can use 
the double coating to produce a manifold with a different topology, but with 
the change in topology invisible to external measurements. 

To define the solutions of Maxwell's equations rigorously in the single coating 
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case, one has to add boundary conditions on £. Physically, this corresponds 
to the lining of the interior of the single coating material, e.g., in the case 
of blowing up a point, with a perfectly conducting layer, see j3]. We point 
out here that in the recent preprint [21], the single coating construction is 
supplemented with selfadjoint extensions of Maxwell operators in the interior 
of the cloaked regions; these implicitly impose interior boundary conditions 
on the boundary of the cloaked region, similar to the PEC boundary condi- 
tion suggested in [3]. For the case of an infinite cylinder the Soft-and-Hard 
(SH) interior boundary condition is used in [3] to guarantee cloaking of active 
objects, and is needed even for passive ones. 

In this paper, we show how more elaborate geometric constructions, cor- 
responding to blowing up a curve, enable the description of tunnels which 
allow the passage of waves between distant points, while only the ends of 
the tunnels are visible to external observation. These devices function as 
electromagnetic wormholes, essentially changing the topology of space with 
respect to solutions of Maxwell's equations. 

We form the wormhole device around an obstacle K C M 3 as follows. First, 
one surrounds K with metamaterials, corresponding to a specification of EM 
parameters e and fx. Secondly, one lines the surface of K with material im- 
plementing the Soft-and-Hard (SH) boundary condition from antenna theory 
[3 HOI EH]; this condition arose previously [3] in the context of cloaking an 
infinite cylinder. The EM parameters, which become singular as one ap- 
proaches K, are given as the pushforwards of nonsingular parameters e and 
fi on an abstract three-manifold M, described in Sec. 2. For a curve 7 C M, 
we construct the diffeomorphism F from M\j to the wormhole device in Sec. 
3. For the resulting EM parameters e and fl, we have singular coefficients of 
Maxwell's equations at K, and so it is necessary to formulate an appropriate 
notion of locally finite energy solutions (see Def. 4.1). In Theorem 4.2, we 
then show that there is a perfect correspondence between the external mea- 
surements of EM waves propagating through the wormhole device and those 
propagating on the wormhole manifold. 

It was shown in [3] that the cloaking constructions are mathematically valid 
at all frequencies k. However, both cloaking and the wormhole effect stud- 
ied here should be considered as essentially monochromatic, or at least very 
narrow-band, using current technology, since, from a practical point of view 
the metamaterials needed to implement the constructions have to be fabri- 
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cated and assembled with a particular wavelength in mind, and theoretically 
are subject to significant dispersion [TS]. Thus, as for cloaking in [Tlfl [TSl 13] . 
here we describe the wormhole construction relative to electromagnetic waves 
at a fixed positive frequency k. We point out that the metamaterials used in 
the experimental verification of cloaking [20] should be readily adaptable to 
yield a physical implementation, at microwave frequencies, of the wormhole 
device described here. See Remark 1 in Sec. 14.21 for further discussion. 

The results proved here were announced in [3]. 

2 The wormhole manifold M 

First we explain, somewhat informally, what we mean by a wormhole. The 
concept of a wormhole is familiar from general relativity [HJ [22], but here 
we define a wormhole as an object obtained by gluing together pieces of 
Euclidian space equipped with certain anisotropic EM parameter fields. We 
start by describing this process heuristically; later, we explain more precisely 
how this can be effectively realized vis-a-vis EM wave propagation using 
metamaterials. 

We first describe the wormhole as an abstract manifold M, see Fig. 1; in 
the next section we will show how to realize this concretely in M 3 , as a 
wormhole device N. Start by making two holes in the Euclidian space M 3 = 
{(x,y,z)\x,y,z e E}, say by removing the open ball B\ = B(0,1) with 
center at the origin O and of radius 1, and also the open ball B 2 = B(P, 1), 
where P = (0, 0, L) is a point on the z-axis having the distance L > 3 to 
the origin. We denote by Mi the region so obtained, Mi = IR 3 \ (B\ U B 2 ), 
which is the first component we need to construct a wormhole. Note that 
Mi is a 3-dimensional manifold with boundary, the boundary of Mi being 
dMi = dBi U dB 2 , the union of two 2-spheres. Thus, dM% can be considered 
as a disjoint union § 2 US 2 , where we will use S> 2 to denote various copies of 
the two-dimensional unit sphere. 

The second component needed is a 3— dimensional cylinder, M 2 = § 2 x [0, 1]. 
This cylinder can be constructed by taking the closed unit cube [0, l] 3 in M 3 
and, for each value of < s < 1, gluing together, i.e., identifying, all of the 
points on the boundary of the cube with z = s. Note that we do not identify 
points at the top of the boundary, at z = 1, or at the bottom, at z = 0. We 



4 



then glue together the boundary dB(0, 1) ~ § 2 of the ball B(0, 1) with the 
lower end (boundary component) S> 2 x {0} of M 2 , and the boundary dB(P, 1) 
with the upper end, § 2 x {1}. In doing so we identify the point (0,0, 1) G 
dB(0,l) with the point NP x {0} and the point (0, 0, L - 1) G OB (P, 1) 
with the point NP x {1}, where NP is the north pole on §> 2 . 

The resulting manifold M no longer lies in R 3 , but rather is the connected 
sum of the components M\ and M 2 and has the topology of M 3 with a 
3— dimensional handle attached. Note that adding this handle makes it pos- 
sible to travel from one point in M\ to another point in Mi, not only along 
curves lying in M\ but also those in M 2 . 

To consider Maxwell's equations on M, let us start with Maxwell's equations 
on M. 3 at frequency k G M., given by 

VxE = ikB, VxH = -ikD, D(x) = e(x)E(x), B(x) = fi(x)H(x). 

Here E and H are the electric and magnetic fields, D and B are the electric 
displacement field and the magnetic flux density, e and /i are matrices corre- 
sponding to permittivity and permeability. As the wormhole is topologically 
different from the Euclidian space 1R 3 , we use a formulation of Maxwell's 
equations on a manifold, and as in [3], do this in the setting of a general Rie- 
mannian manifold, (M,g). For our purposes, as in [TU [3] it suffices to use 
e, fi which are conformal, i.e., proportional by scalar fields, to the metric g. 
In this case, Maxwell's equations can be written, in the coordinate invariant 
form, as 

dE = ikB, dH = -ikD, D = eE, B = fiH in M, 

where E, H are 1-forms, D, B are 2-forms, d is the exterior derivative, and e 
and fi are scalar functions times the Hodge operator of (M,g), which maps 
1-forms to the corresponding 2-forms [2]. In local coordinates these equations 
are written in the same form as Maxwell's equations in Euclidian space with 
matrix valued e and /i. Although not necessary, for simplicity one can choose 
a metric on the wormhole manifold M which is Euclidian on Mi , and on M 2 
is the product of a given metric go on S 2 and the standard metric of [0, 1]. 
More generally, can also choose the metric on M 2 to be a warped product. 
Even the simple choice of the product of the standard metric of S 2 and the 
metric 5 2 ds 2 , where 5 is the "length" of the wormhole, gives rise to interesting 
ray-tracing effects for rays passing through the wormhole tunnel. For 5 << 1, 
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the image through one end of the wormhole (of the region beyond the other 
end) would resemble the image in a a fisheye lens; for 8 > 1, multiple images 
and greater distortion occur. (See [U Fig.2].) 

The proof of the wormhole effect that we actually give is for yet another 
variation, where the balls that form the ends have their boundary spheres 
flattened; this may be useful for applications, since it allows for there to be a 
vacuum (or air) in a neighborhood of the axis of the wormhole, so that, e.g., 
instruments may be passed through the wormhole. We next show how to 
construct, using metamaterials, a device N in R 3 that effectively realizes the 
geometry and topology of M, relative to solutions of Maxwell's equations at 
frequency k, and hence functions as an electromagnetic wormhole. 

3 The wormhole device N in IR 3 

We now explain how to construct a "device" N in M 3 , i.e., a specification 
of permittivity e and permeability fi, which affects the propagation of elec- 
tromagnetic waves in the same way as the presence of the handle M 2 in the 
wormhole manifold M. What this means is that we prescribe a configuration 
of metamaterials which make the waves behave as if there were an invisible 
tube attached to IR 3 , analogous to the handle M 2 in the wormhole manifold 
M. In the other words, as far as external EM observations of the wormhole 
device are concerned, it appears as if the topology of space has been changed. 

We use cylindrical coordinates (9, r, z) corresponding to a point (r cos 9, r sin 9, z) 
in IR 3 . The wormhole device is built around an obstacle K C R 3 . To de- 
fine K, let 5* be the two-dimensional finite cylinder {9 G [0,27r],r = 2, < 
z < L} C IR 3 . The open region K consists of all points in IR 3 that have 
distance less than one to S and has the shape of a long, thick-walled tube 
with smoothed corners. 

Let us first introduce a deformation map F from M to N = IR 3 \ K or, more 
precisely, from M\j to iV\£, where 7 is a closed curve in M to be described 
shortly and S = OK. We will define F separately on M\ and M 2 denoting 
the corresponding parts by F\ and F 2 . 

To describe F\, let 71 be the line segment on the z— axis connecting dB(0, 1) 
and dB(P, 1) in M u namely, 71 = {r = 0, z G [1, L - 1]}. Let F x (r,z) = 
(9, R(r, z), Z(r, z)) be such that (i?(r, z), Z(r, z)), shown in Fig. 2, 
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Figure 1: Schematic figure: a wormhole manifold is glued from two com- 
ponents, the "handle" and space with two holes. Note that in the actual 
construction, the components are three dimensional. 




~A B C D A' D' 



Figure 2: The map (R(r, z) , Z (r, z)) in cylindrical coordinates (z,r). 
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transforms in the (r, z) coordinates the semicircles AB and CD m. the left 
picture to the vertical line segments A'B' — {r G [0, 1], z — 0} and CD' = 
{r G [0, 1], z = L] in the right picture and the cut 71 on the left picture to the 
curve B'C on the right picture. This gives us a map F\ : M\ \ 71 — ► N\ \ E, 
where the closed region iVi in R 3 is obtained by rotation of the region exterior 
to the curve A'B' CD' around the z— axis. We can choose F\ so that it is 
the identity map in the domain [/ = M 3 \ { — 2 < z < L + 2, < r < 4}. 

To describe F 2 , consider the line segment, 72 = {NP} x [0, 1] on M 2 . The 
sphere without the north pole can be "flattened" and stretched to an open 
disc with radius one which, together with stretching [0, 1] to [0, L], gives us a 
map F 2 from M 2 \7 2 to iV 2 \E. The region N 2 is the 3— dimensional cylinder, 
N 2 = {6 G [0,2vr],r G [0,1], z G [0,L]}. When flattening § 2 \ NP, we do 
it in such a way that F\ on dB(0,l) and dB(P,l) coincides with F 2 on 
(§ 2 \ NP) x {0} and (§ 2 \ NP) x {1}, respectively. 

Thus, F maps M \ 7, where 7 = 71 U 7 2 is a closed curve in M, onto iV \ E; 
in addition, F is the identity on the region U. 

Now we are ready to define the electromagnetic material parameter tensors 
on N. We define the permittivity to be 



Fjf£{y)= iDF)(x)-e(x)-(DF(x)y 



det(DF) 



--F-Hy) ' 



where DF is the derivative matrix of F, and similarly the permeability to 
be /x = F*/jl. These deformation rules are based on the fact that permittivity 
and permeability are conductivity type tensors, see 



Maxwell's equations are invariant under smooth changes of coordinates. This 
means that, by the chain rule, any solution to Maxwell's equations in M \ 7, 
endowed with material parameters e, // becomes, after transformation by F, 
a solution to Maxwell's equations in N \ E with material parameters e and ju, 
and vice versa. However, when considering the fields on the entire spaces M 
and N, these observations are not enough, due to the singularities of e and Jl 
near E; the significance of this for cloaking was observed and analyzed in [3]. 
In the following, we will show that the physically relevant class of solutions to 
Maxwell's equations, namely the (locally) finite energy solutions, remains the 
same, with respect to the transformation F, in (M;e,ji) and (N;e,Jl). One 
can analyze the rays in M and N endowed with the electromagnetic wave 
propagation metrics g = y/ejl and g = respectively. Then the rays on 
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M are transformed by F into the rays in N. As almost all the rays on M do 
not intersect with 7, therefore, almost all the rays on N do not approach E. 
This was the basis for [T6j [T8] and was analyzed further in [19] ; see also [T7] 
for a similar analysis in the context of elasticity. Thus, heuristically one is 
led to conclude that the electromagnetic waves on (M; e, ll) do not feel the 
presence of 7, while those on (AT; e, /I) do not feel the presence of K, and 
these waves can be transformed into each other by the map F. 

Although the above considerations are mathematically rigorous, on the level 
both of the chain rule and of high frequency limits, i.e., ray tracing, in the 
exteriors M \ 7 and iV \ S, they do not suffice to fully describe the behavior 
of physically meaningful solution fields on M and N. However, by carefully 
examining the class of the finite-energy waves in M and N and analyzing 
their behavior near 7 and S, respectively, we can give a complete analysis, 
justifying the conclusions above. Let us briefly explain the main steps of the 
analysis using methods developed for theory of invisibility (or cloaking) at 
frequency k > [3] and at frequency k = in [U [7] . The details will follow. 

First, to guarantee that the fields in N are finite energy solutions and do not 
blow up near S, we have to impose at £ the appropriate boundary condition, 
namely, the Soft-and-Hard (SH) condition, see [SHU], 

e e -£| s = 0, eg-H\ E = 0, 

where eg is the angular direction. Secondly, the map F can be considered 
as a smooth coordinate transformation on M \ 7; thus, the finite energy 
solutions on M \ 7 transform under F into the finite energy solutions on 
N \ E, and vice versa. Thirdly, the curve 7 in M has Hausdorff dimension 
equal to one. This implies that the possible singularities of the finite energy 
electromagnetic fields near 7 are removable J12J, that is, the finite energy 
fields in M \ 7 are exactly the restriction to M \ 7 of the fields defined on all 
of M. 

Combining these steps we can see that measurements of the electromagnetic 
fields on (M; e, /1) and on {R 3 \ K; e, Jm) coincide in U. In the other words, if 
we apply any current on U and measure the radiating electromagnetic fields 
it generates, then the fields on U in the wormhole manifold (M; e, //) coincide 
with the fields on U in (M 3 \ K; e, fx) , 3-dimensional space equipped with the 
wormhole device construction. 

Summarizing our construction, the wormhole device consists of the metama- 
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terial coating of the obstacle K. This coating should have the permittivity 
e and permeability Jl. In addition, we need to impose the SH boundary 
condition on E, which may be realized by fabricating the obstacle K from a 
perfectly conducting material with parallel corrugations on its surface [H [11] . 

In the next section, the permittivity e and and permeability Jl are described 
in a rather simple form. (As mentioned earlier, in order to allow for a tube 
around the axis of the wormhole to be a vacuum or air, we deal with a 
slightly different construction than was described above, starting with flat- 
tened spheres). It should be possible to physically implement an approxima- 
tion to this mathematical idealization of the material parameters needed for 
the wormhole device, using concentric rings of split ring resonators as in the 
experimental verification of cloaking obtained in [20J. 

4 Rigorous construction of the wormhole 

Here we present a rigorous model of a typical wormhole device and justify 
the claims above concerning the behavior of the electromagnetic fields in the 
wormhole device in R 3 in terms of as the fields on the wormhole manifold 
(M,g). 

4.1 The wormhole manifold (M,g) and the wormhole 
device N 

Here we prove the wormhole effect for a variant of the wormhole device 
described in the previous sections. Instead of using a round sphere S 2 as 
before, we present a construction that uses a deformed sphere §g at that is 
flat the near the south and north poles, SP and NP. This makes it possible 
to have constant isotropic material parameters near the z-axis located inside 
the wormhole. For possible applications, see [4]. 

We use following notations. Let (9, r, z) G [0, 2n] x R + x R be the cylindrical 
coordinates of R 3 , that is the map 

X : (9, r, z) — > (r cos 9, r sin 9, z) 

that maps X : [0, 2vr] x R + x R -> R 3 . In the following, we identify [0, 2tt] 
and the unit circle S 1 . 
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Let us start by removing from M 3 two "deformed" balls which have flat 
portions near the south and north poles. More precisely, let Mi = M 3 \ (Pi U 
P2), where in the cylindrical coordinates 

P 1 = {X(9, r,z) : -1 < z < 1, < r < 1} 

Li{X(9,r,z) : (r- if + z 2 < 1}, 
P 2 = {X(6,r,z) : 1 < z-L < 1, < r < 1} 

U{X(0,r,z): (r-l) 2 + (^-L) 2 < 1}. 

We say that the boundary <9Pi of Pi is a deformed sphere with flat portions, 
and denote it by Sg at . We say that the intersection points of §g at with the 
z-axis are the north pole, NP, and the south pole, SP. 

Let gi be the metric on Mi inherited from M 3 , and let 71 be the path 
7l = {X(0, 0, z) : 1 < z < L - 1} C Mi. 

Set 

Ai = Mi \ Vi /4 , 

VJ = {X(^,r,^) : < r < t, 1 < z < L- 1}, < t < 1, 

and consider a map G : Mi \ f y 1 — > see Fig. 3. G defined as the identity 
map on Mi \ V1/2 and, in cylindrical coordinates, as 

G (X(6, r, z)) = X(6, \ + (6,r,z)e V 1/2 . 

Clearly, G is C ' 1 — smooth. 

Let U(x) G R 3x3 , x = X(6,r,z), be the orthogonal matrix that maps the 
standard unit vectors ei, e 2 , e 3 of 1R 3 to the Euclidian unit vectors correspond- 
ing to the 6, r, and z directions, that is, 

U(x)ei = (- sin 0, cos 0,0), U(x)e 2 = (cos 6, sin 6, 0), C/(x)e 3 = (0, 0, 1). 

Then the differential of G in the Euclidian coordinates at the point x £ V1/2 
is the matrix 

DG (a:)f/(y) § CZ(x)- 1 , x = X(0,r,z), y = G (x). (1) 

V 1 y 
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Later we impose on part of the boundary, S = dA\ n{l<z<L — 1}, the 
soft-and-hard boundary condition (marked red in the figures). 

Next, let (9,z,r) = (0(x), z(x), t{x)) be the Euclidian boundary normal 
coordinates associated to So, that is, t(x) = distR3(x, S ) and (9(x), z(x)) 
are the 9 and ^-coordinates of the closest point of So to x. 

Denote by (Cro)*<7i the push forward of the metric g\ in Go, that is, the 
metric obtained from gi using the change of coordinates Go, see [2J. The 
metric (Go)*gi coincides with g x in A x \ V\i%, and in the Euclidian boundary 
normal coordinates of S , on A\ fl Vj/2, the metric (Go)*<7i, has the length 
element 

ds 2 = At 2 d6 2 + dz 2 + Adr 2 . 





Go 







' ) 




1 





Figure 3: A schematic figure on the map Go, considered in the (r, z) coor- 
dinates. Later, we impose the SH boundary condition on the portion of the 
boundary coloured red. 

Next, let 

<?3 = conv({(r,^): (r - 2) 2 + (z - (-2)) 2 < 1} 

U{(r,z): (r-2) 2 + (z-(L + 2)) 2 <l}), 
q A = {( r ,z) : < r < 1, -1 < z < L+ 1}, 
where conv(g) denotes the convex hull of the set q. 
Let 

Ni = M 3 \(P 3 UP 4 ), 
P 3 = {X{6,r,z): (r,z)eq 3 }, 
P A = {X{6,r,z): (r,z)eq,}, 
S x = dNiXdPi. 
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We can find a Lipschitz smooth map Gx : Ax — > N ± , see Fig. 4, of the form 

G 1 (X(e,r,z))=X(6,R(r,z),Z(r,z)) 

such that it maps S to Si, and in A 1 near S it is given by 

G 1 (x + tu ) = G 1 (x)+tu 1 . (2) 

Here, x G S , vq is the Euclidian unit normal vector of S , v\ is the Euclidian 
unit normal vector of Si, and < t < \. Moreover, we can find a Gx so that 
it is the identity map near the z-axis, that is, 

G x {x) = x, x G A x n {0 < r < i} (3) 

and such that Gx is also the identity map in the set of points with the 
Euclidian distance 4 or more from Pi U P<i. Note that we can find such a 
Gx such that both Gx and its inverse G^ 1 are Lipschitz smooth up to the 
boundary. Thus the differential DGx of Gx at x G Ax in Euclidian coordinates 
is 

DGx(x) = U(y)(^ ail ^ z) A ^ z) )u{x)-\ x = X(6, r, z), y = G 1 (x), 

where cq < an(r, z) < ci and A(r, z) is a symmetric (2 x 2)-matrix satisfying 

col < A(r, z) < cxl 

with some co,cx > 0. 

The map Fx(x) = Gx(G (x)) then maps Fx : Mx\"/x — > iVi. Let #i = (Fx)*gx 
be metric on JVi. From the above considerations, we see that the differential 
DFx of Fx at x G M x \ 71 near S , in Euclidian coordinates, is given by 

DF 1 ( x )=U(y)^' > - r ' z) B ^ z) )u(x)-\ (4) 

W-'-'^ mJ™?,*),*) x = x ^ » = 

where c < Cn(r, 2;) < Ci, and B(r, z) is a symmetric (2 x 2)-matrix satisfying 
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for some c , c\ > 0. 

Note that <9P 4 n{r < 1} consists of two two-dimensional discs, -82(0, 1) x { — 1} 
and -62(0, 1) x {L + 1}. Below, we will use the map 

h = Fi\ dPl \ N p : dP x \ NP -> B 2 (0, 1) x {-1} C dN x . 

The map f 2 can be considered as the deformation that "flattens" Sg at \ NP 
to a two dimensional unit disc. 




Figure 4: Map G\ in (r, ^-coordinates. 

To describe f 2 , consider §^ at as a surface in Euclidian space and define on it 
the 9 coordinate corresponding to the 9 coordinate of M 3 \ {z = 0}. Let then 
s(y) be the intrinsic distance of y G Sg at to the south pole SP. Then (9, s) 
define coordinates in Sg Bt \{S'P, A^P}. We denote by y(0, s) G S^ at \{^P, A^P} 
the point corresponding to the coordinates (9, s). 

By the above construction, the map f 2 has the form, with respect to the 
coordinates used above, 

f 2 (y(9,s))=X(9,R(s),-l)eB 2 (0,l)x{-l}, where (5) 

R(s) = s, for < s < 

P(s) = 1 - -[(tt + 4) - s], for (vr + 4) - i < s < (vr + 4), 

cf. formulae (j2J) and ([3]). In the following we identify B 2 (0, 1) x { — 1} with 
the disc B 2 (Q, 1). 

Let /ii be the metric on dP\ \ NP inherited from (M\,gx). Let h 2 = (f 2 )*hi 
be the metric on B 2 (0,1). We observe that the metric h 2 makes the disc 
P2(0, 1) isometric to Sg at \ NP, endowed with the metric inherited from K 3 . 
Thus, let 

M 2 =§L t x[-l,I+l]. 
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On M 2 , let the metric g 2 be the product of the metric of §g at inherited from IR 3 
and the metric a 2 (z)dz 2 , a 2 > on [— 1, L + 1]. Let 72 = {-/VP} x [— 1, L + l] 
be a path on M 2 . 

Define jV 2 = P 4 = {X(0, r, *) : < r < 1, -1 < z < L + 1} C M 3 , 
S 2 = 9A^2 H {r = 1}, and let F 2 : M 2 \ 72 — > N 2 be the map of the form 

F 2 {y, z) = {f 2 (y), z) e M 3 , (y, z) E (§L \ NP) x [—1, L+l]. (6) 

Let g 2 = (F 2 )*g 2 be the resulting metric on A^. 



\ 1 

Figure 5: The set N 2 in the (r, z) coordinates. Later, we impose the SH 
boundary condition on the portion of the boundary colored red. 

Denote by M 1 = MiU<9Mi_the closure of M x and let (M, g) = (M 1 ,g 1 )^(M 2 , g 2 ) 
be the connected sum of Mi and M 2 , that is, we glue the boundaries dMi 
and dM 2 . The set N = N% U N 2 C M 3 is open, and its boundary <9iV is 
S = Ei U E 2 - 

Let F be the map F : M \ 7 — > N defined by the maps F\ : Mi \ 71 — > Nx 
and F 2 : M 2 \ 72 — > N 2 , and finally, let 7 = 71 U 72 and g = F*g. 




Figure 6: The set N = Ni U N 2 C M 3 having the complement i^, presented 
in the (r, z) coordinates. Later, the SH boundary condition is imposed on 
OK. 

Let K = M. 3 \ N. On the surface £ = OK we can use local coordinates (t, 6), 
where 9 is the ^-coordinate of the ambient space M 3 and t is either the r or 
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z -coordinate of the ambient space M. 3 restricted to S. Denote also 

r = t(x) = dist]g3(x, dK). 

Then by formula ([2]) we see that in Ni, in the Euclidian boundary normal 
coordinates (9,t,r) associated to the surface Si, the metric g has the length 
element 

ds 2 = Adr 2 + dt 2 + At 2 d9 2 , < r < ~, c 1 < a^t) < c , c > 1. 

The construction of _F 2 yields that in N 2 , in the Euclidian boundary normal 
coordinates (9,t,r) with t — z, associated to the surface S 2 = dK fl <9iV2, 
the metric g has the length element, near E 2 , 

ds 2 = 4d^ + a 2 (t)dP + 4¥ 2 d9 2 , 0<r<^. 

Here, near 9A^i fl <9iV2, we use t = z on Si. Choosing the map G\ in the 
construction of the map F x appropriately, we have a 2 (— 1) = cni(— 1), a 2 (L + 
1) = «i(L + 1), and the resulting map is Lipschitz. 

On Mi, Ni, and N 2 that are subsets of M. 3 we have the well defined cylin- 
drical coordinates (9, r, z). Similarly, M 2 = §g at x [— 1,L + 1] we define 
the coordinates (9, s, z), where (9, s) are the above defined coordinates on 
§l t \{SP,NP}. 

We can also consider on iV C M. 3 also the Euclidian metric, denoted by 
g e . In Euclidean coordinates, (g e )ij = Sjk- Consider next the above defined 
Euclidian boundary normal coordinates (9, t, r) associated to dK. They are 
well defined in a neighborhood of dK. We define the vector fields 

on N near dK. These vector fields are orthogonal with respect to the metric 
g and to the metric g e . 

On M near 7, we use coordinates (9, t, r). On Mi, near 71 they in the terms 
of the cylindrical coordinates are (9,t,r) = (9,z,r). On M 2 , they are the 
coordinates (9,t,r) = (9,z,s), where s is the intrinsic distance to the north 
pole NP. We define also the vector fields 

£ = d T , r) = d g , C = d t 
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on M \ 7 near 7. These vector fields are orthogonal with respect to the 
metric g. 

In the sequel, we consider the differential of F as the linear map 

DF : (T x M } g) - (T y N,g e ), y = F(x), x G M \ 7. 

Using formula (j3J) in Mi and formulas (jSJ), ([6]) in M2, we see that .DF -1 ^) 
at x € iV near dN is a bounded linear map that satisfies 

(7) 

\^DF-\x)?f) g \ <Cr(x) } (C,DF-\x)rj) g = 0, (£, DF^x)^ = 0, 
{^DF-\x)C) a = 0, KCDF-^xXU < C, \(Z,DF-\x)C) g \ < C, 
(77, DF" 1 ^ = 0, KC^H^ai < C, K^^^OsI < C, 

where C > and (• , • ) s is the inner product defined by the metric g. More- 
over, we obtain similar estimates for DF in terms of the Euclidian metric 

9\ 



\(%DF(y)r ] ) g ,\ < Cr(y)-\ ((, DF(y) V ) g e = 0, (^DF(y) V ) g e = 0, 
{%DF{y)Q g e = 0, \(C,DF{y)CU < C, \(£,DF(y)() g e\ < C, 
(%DF(y)§ ge = 0, \(C,DF(y)£) ge \ < C, \(tDF(y)£) ge \ < C 

for y G M \ 7 near 7 with C > 0. 

Next, consider DF(y) at y G M\j. Recall that the singular values Sj(y), j = 
1,2,3 of DF(y) are the square roots of the eigenvalues of (DF(y)) t DF(y), 
where {DF) 1 is the transpose of DF. By ([7]), the singular values Sj = Sj(y), 
j = 1, 2, 3, of DF(y), numbered in increasing order, satisfy 

ci < si(y) < c 2 , 
ci < s 2 (y) < c 2 , 

< 8 3 {y) < 



r(y) r(yY 
where Ci, c 2 > 0. 

The determinant of the matrix DF(y) can be computed in terms of its sin- 
gular values by det (DF) = Sis 2 s 3 . Later, we need the norm of the matrix 
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det(DFiy))- 1 DF(y). It satisfies by formula © 

3 

WdetiDFiy))- 1 DF(y)\\ = IKjJ^diag (s u s 2} s 3 )\\ = max J]^ 1 < cf. (9) 

k=l 3 k^j 

4.2 Maxwell's equations on the wormhole with SH coat- 
ing 

Let dVo(x) denote the Euclidian volume element on iV C 1R 3 . Recall that 
Af C I 3 is open set with boundary dN = E. Let dV g be the Riemannian 
volume on (M, g). We consider below the map F : M\^ — > N as a coordinate 
deformation. The map F induces for any differential form E on N a form 
E = F*E in M \ 7 called the pull back of E in F, see [2]. 

Next, we consider Maxwell equations with degenerate material parameters e 
and fx on N with SH boundary conditions on E. On M and iV we define the 
permittivity and permeability by setting 

£ jk = ^jk = djet ( gr /2gik i on M , (10) 

e> k = J,i k = det(g) 1/2 g jk , on N. 

Here, and below, the matrix [gjk{x)} is the representation of the metric g in 
local coordinates, [g^ k {x)] is the inverse of the matrix [gjk(x)}, and det(g) is 
the determinant of [gjk(x)}- We note that the metric g is degenerate near 
E, and thus e and jl, represented as matrices in the Euclidian coordinates, 
have elements that tend to infinity at E, that is, the matrices e and ]1 have 
a singularity near E. 

Remark 1. Modifying the above construction by replacing M 2 with M 2 = 
§jj at x l 2 ] for appropriate fa, l 2 G M and choosing in an appropriate way, 
we can use local coordinates (6, t) on E such that the Euclidian distance 
along E of points (6, fa) and (9,t 2 ) is proportional to \ fa — t 2 \, and the metric 
g in the Euclidian boundary normal coordinates (9, t, r) associated to dK 
has the form 

ds 2 = 4rfr 2 + dt 2 + 4r 2 d9 2 , < r < -. 

4 

The metric corresponding to the metamaterials used in the physical exper- 
iment in [20] has the same form in Euclidian boundary normal coordinates 
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associated to an infinitely long cylinder -82(0, 1) x !L Thus it seems likely that 
metamaterials similar to those used in the experimental verification of cloak- 
ing could be used to create physical wormhole devices working at microwave 
frequencies. 



4.3 Finite energy solutions of Maxwell's equations 
and the equivalence theorem 

In the following, we consider 1-forms E = £\ Ejdx j and H = J2j Hjdx j in 
the Euclidian coordinates (x 1 , rc 2 , J 3 ) of N C R 3 . In the sequel, we use Ein- 
stein's summation convention and omit the sum signs. We use the Euclidian 
coordinates as we want to consider N with the differential structure inherited 
from the Euclidian space. We say that Ej and Hj are the (Euclidian) coeffi- 
cients of the forms E and H, correspondingly. We say that these coefficients 
are in L p oc (N ', dV ), 1 < p < 00, if 



/ 

Jw 



\Ej(x)\ p dV (x) < 00, for all bounded measurable sets W C N. 



Definition 4.1 We say that the 1-forms E and H are finite energy solutions 
of Maxwell's equations in N with the soft-and-hard (SH) boundary conditions 
on S and the frequency k ^ 0, 

Vx£ = ikJl(x)H, V x H = -ike(x)E + J on N, 
77-£| E = 0, ?j-H\x = 0, 

if 1-forms E and H and 2-forms D = eE and B = JiH in N have coefficients 
in Ll oc (N,dV ) and satisfy 

\\E\\l*{w,\g\ 1 '*W>)= / ^ jk E j E k dV (x)<oo, 
11^11^(^,1/2^)) = J JP k HjH k dV {x) < 00 
for all bounded measurable sets W C N, and finally, 

[ ((V x h) ■ E - ikh ■ Jl(x)H) dV (x) = 0, 

JN 

[ ((V x e) • H + e ■ {ike{x)E - J)) dV (x) = 0, 

JN 



19 



for all 1-forms e and h with coefficients in C£°(N) that satisfy 

V-e\x = 0, rj-h\x = 0, (11) 
where rj = dg is the angular vector field that is tangential to E. 

Below, we use for 1-forms E = Ejdx^ and H = Hjdx^ , given in local coordi- 
nates (x 1 , x 2 , x 3 ) on M, the notations 

Vx£ = dH, V- (eE) = d*E, V- (fiH) = d*H, 

where d is the exterior derivative and * is the Hodge operator on (M,g), cf. 
formula ffTUj) . 

We have the following "equivalent behavior of electromagnetic fields on N 
and M" result, analogous to the results of [3] for cloaking. 

Theorem 4.2 Let E and H be 1-forms on M \ 7 and E and H be 1-forms 
with coefficients in Ll oc (N,dV ) such that E = F*E, H = F*H. Let J 
and J = F*J be 2- forms with smooth coefficients in N and M \ 7 that are 
supported away from E and 7 . 

Then the following are equivalent: 

1. On N , the 1-forms E and H satisfy Maxwell's equations with SH bound- 
ary conditions in the sense of Definition al 

2. On M , the forms E and H can be extended on M so that they are 
classical solutions E and H of Maxwell's equations, 

V x E = ikfiH, in M, 

V x H = -ikeE + J, %n M. 

Proof. Assume first that E and H satisfy Maxwell's equations on M with 
source J supported away from 7. Then E and H are C°° smooth near 7. 

Using F _1 : iV — > M \ 7 we define the 1-forms E, H and 2-form J on N 
by E = (F- X )*E, H = (F- X )*H, and J = (F -1 )* J. These fields satisfy 
Maxwell's equations in N, 

V x E = ikJl(x)H, V x H = -ike(x)E + J in N. (12) 
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Now, writing E = Ej{x)dx^ on M near 7, we see using the transformation 
rule for differential 1-forms that the form E = (F^yE is in local coordinates 



E = Ej(x)dx j , Ej(x) = {DF~ l )){x) E k (F'\x)), x e N. (13) 



Using the smoothness of E and H near 7 on M and formulae (j7|), we see 
that E, H are forms on N with L\ oc {N, dVo) coefficients. Moreover, 



where x G N, y = F G M \ 7. Formula (Q shows that D = eE, and 
B = JiH are 2- forms on iV with L\ oc (N, dVo) coefficients. 

Let C N be the i-neighbourhood of E in the (/-metric. Note that for 
small t > the set E(t) is the Euclidian (t/2)-neighborhood of OK. Denote 
by v be the unit exterior Euclidian normal vector of <9E(t) and the Euclidian 
inner product by (rj, E) g e —rj-E. 

Formulas ([7j) and fflBl imply that the angular components satisfy 



with some C > 0. Thus denoting by dS the Euclidian surface area on <9E(t), 
Stokes' formula, formula ffl2|) . and the identity v x ^ = ±rj yield 



det(DF(y))- 1 J DF( 2 /)£(y) J DF( 2/ ) t ( J DF( 2/ ) t )- 1 J E(y) 
det(DF(y)r 1 DF(y)e(y)E(y) 



\r)-E\<Ct, xedZ(t) 



and 



|C--B|<C, xG<9E(t) 
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for a test function h satisfying formula (II ip . 

Similar analysis for H shows that 1-forms E and H satisfy Maxwell's equa- 
tions with SH boundary conditions in the sense of Definition 14.11 

Next, assume that E and H form a finite energy solution of Maxwell's equa- 
tions on (N, g) with a source J supported away from E, implying in particular 
that 

P k E-E k eL\W,dVo), p k H~H~ k EL\W, dV ) 

where W = F(U \ 7) C iV and U C M is a relatively compact open neigh- 
bourhood of 7, supp (J)nW = 0. Define E = F*E, H = F*H, and J = F*J 
on M \ 7. Therefore we conclude that 

V x E = ik/i(x)H, V x H = -ike{x)E + J, in M \ 7 

and 

e jk E^W h eL\U\ 7, dV g ), p? k E^h G L\U \ 7, dV g ). 

As representations of e and /x, in local coordinates of M, are matrices that 
are bounded from above and below, these imply that 

VxEe L 2 (U \ 7, dV g ), VxH G L 2 (U \ 7, dV g ), 
V- (eE) = 0, V- (jiH) = 0, in £/ \ 7. 

Let .E 6 , if e G L 2 (U, dVg) be measurable extensions of E and if to 7. Then 

V x £ e - ikfi(x)H e = 0, fnC/\7, 

V x E e - tkfi(x)H e G if -1 ^, dV g ), 

V x H e + ike(x)E e = 0, in f/ \ 7, 
Vx H e + iks(x)E e G if -1 ^, dV^), 

where if _1 (?7, dV^) is the Sobolev space with smoothness (—1) on (U,g). 
Since 7 is a subset with (Hausdorff) dimension 1 of the 3-dimensional domain 
U, it has zero capacitance. Thus, the Lipschitz functions on U that vanish on 
7 are dense in see |12| . Therefore, there are no non-zero distributions 

in H~ X (U) supported on 7. Hence we see that 

VxF- ik/i(x)H e = 0, V x H e + ike(x)E e = in U. 



22 



This also implies that 

V- (eE e ) = 0, V- (/i# e ) = in U, 

which, by elliptic regularity, imply that E e and H e are C°° smooth in U. 

In summary, E and H have unique continuous extensions to 7, and the 
extensions are classical solutions to Maxwell's equations. □ 
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